Breakdown of the Narrow Width Approximation for New Physics 
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The narrow width approximation is used in high energy physics to reduce the complexity of 
scattering calculations. It is a fortunate accident that it works so well for the Standard Model, 
but in general it will fail in the context of new physics. We find numerous examples of significant 
corrections when the calculation is performed fully off-shell including a finite width, notably from 
effects from the decay matrix elements, not just phase space. If not taken into account, attempts 
to reconstruct the Lagrangian of a new physics discovery from data would result in considerable 
inaccuracies and likely inconsistencies. 



The narrow-width approximation (NWA) is used ex- 
tensively to calculate cross sections for production of 
promptly-decaying particles. It's use is justified only if 
five critical conditions are met: (i) the total width of a 
resonant particle is much smaller than its mass, T <C M; 
(ii) daughter particles are much less massive than the par- 
ent, m <C M; (hi) the scattering energy is much larger 
than the parent mass, yfs ~3> M; (iv) there is no signif- 
icant interference with non-resonant processes; and (v) 
the resonant propagator is separable from the matrix el- 
ement (ME). If these are valid, the propagator can be 
integrated independently over all q 2 (including unphysi- 
cal values, with negligible effect) to obtain a constant: 
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In a nutshell, the NWA assumes that the massive state is 
always produced exactly at its pole as an asymptotic final 
state, so its decay is an independent process, expressed 
by a simple numerical constant known as a branching ra- 
tio (BR) : the fractional probability to decay to a specific 
final state. Parametrically, the NWA introduces an esti- 
mated error of 0(T/M). The NWA is widely used when 
r/M is small, regardless of the other conditions. While 
there's some awareness of (i), (iii) and (iv), assumptions 
(ii) and (v) have not previously been discussed in the 
literature or textbooks to the best of our knowledge. 

The NWA became standard at a time when numerical 
computation tools were not advanced enough to perform 
off-shell calculations of full matrix elements including de- 
cays and finite widths. Complete analytical calculations 
for such 2 — > n > 2 processes are generally intractable. 
The NWA works well for heavy particle production (W 
and Z bosons and the top quark, excluding hadronic and 
flavor physics) above threshold in the Standard Model 
(SM), largely owing to the fact that the decay products 
are much lighter than the parent. In those cases, T/M 
is indeed small, a couple of percent at most; other un- 
certainties dominate. We need to include finite widths in 
high-energy SM calculations only in a few cases involving 
threshold restrictions, such as scattering at the Z pole, 



e + e~ — * W + W~ — > 4/ production top quark pah- 
production as a background at the Large Hadron Collider 
(LHC) Q, and Higgs boson decay to weak bosons 0]. 
Our work addresses the regime above threshold. 

In the context of new physics extensions to the SM, 
it is generally the case (but not always) that massive 
particle widths are much smaller than their masses, T <C 
M, leading one to conclude that the NWA is still valid. 
Also, collision energy is typically far above production 
threshold, Eqm — M = </s — M » T, thus avoiding a 
cutoff of the Breit-Wigner lineshape. We note a glaring 
exception to this: the proposed technique to measure new 
particle masses via a threshold scan at a future lepton 
collider fij. Only one non-NWA study exists for such a 
case, and for a small subset of particles in one model [f|. 

Rarely does scattering of a given set of initial and final 
states result from only one resonant process. Interference 
with other resonant or non-resonant processes can gener- 
ally occur, rendering the NWA technically inapplicable. 
In the SM at high energy, this is typically insignificant 
compared to other uncertainties, but this is usually not 
true in new physics scenarios. We have found numerous 
instances of significant corrections from interference ef- 
fects. However, these are not the focus of this letter, so 
we defer that discussion to a later work Q . 

Separation of the resonant particle's propagator is, 
strictly speaking, never valid: even when all particles 
are scalars, as the phase space factor for the decay par- 
ticles is a function of q 2 for finite daughter mass, m ^ 0. 
However, for extremely small m the dependence may be 
negligible. The NWA does work very well for top quark 
production and decay, despite a non-trivial matrix ele- 
ment and massive daughter particles. But in general, 
momentum-dependent external wavefunctions and cou- 
plings invalidate Eq. [T] This can result in significant off- 
shell corrections in beyond-the-Standard Model (BSM) 
scenarios, where the NWA is universally adopted. A crit- 
ical ingredient is massive daughters, which most BSM 
scenarios include, often with near-degeneracies driven by 
underlying symmetries. These differences from the SM 
yield unanticipated behavior. 
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Figure 1: Feynman diagrams for various example processes 
discussed in the text: (a) scalar field theory toy model; (b) 
gluino resonance in supersymmetry; (c) VSS and (d) SFF 
momentum-dependent renormalizable interaction vertices. 

We first consider the process of Fig. [T^,, t-channel scat- 
tering in a pure scalar theory. It has no matrix elements 
other than the t- and s-channel propagators. This al- 
lows us to study the s-channel particle's off-shell behav- 
ior without a decay matrix element modification. We 
may later compare these results with those of a specific 
BSM model, where the decay matrix element introduces 
additional momentum-dependent structure. In this toy 
example we have only one possible contributing diagram 
- other topologies may be forbidden by assigning flavor 
or charges to the scalars - and all particles are massless 
except as labeled. The t-channel particle is kept massive 
to avoid a forward-scattering singularity, and only one 
final-state daughter is massive. These simplifications al- 
low us to derive analytical expressions. 

Squaring the matrix element and integrating over 
three-body phase space obtains the 2-^3 off-shell cross 
section. Dividing by the on-shell 2 — ► 2 cross section 
times the naive branching ratio and subtracting 1 ob- 
tains the deviation of the true leading-order rate from 
the NWA result, A = <tofs/&nwa — 1- We defer the 
rather long full result to later work [||. Here we show 
the result for y/s 3> M and finite daughter mass, m^O: 
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where (3 = yl — m 2 /M 2 is the daughter velocity. The 
off-shell result indeed differs from the approximate on- 
shell result by a quantity ostensibly of order T/M (with 
generally mild logarithmic corrections), but with an in- 
verse velocity squared that diverges in the m — » M limit, 
i.e. nearly-degenerate daughter and parent masses. This 
is due to the (j 2 -dependent t-channel propagator and de- 
cay phase space factors, and is our first principal result. 
Its interpretation, however, requires some care. 

The partial width of a scalar particle decaying to two 
scalar daughters is proportional to /? 2 , originating from 



the phase space. If this is the only allowed decay, then the 
total width is equal to the partial width and the velocity 
factors cancel. This leaves only the log terms as correc- 
tion factors to a coefficient which is of order a, x m/M, 
where a t — g 2 /An is the coupling strength of the decay 
interaction. In such a scenario, the leading-order off-shell 
calculation can give 0(cti) rate corrections for decays to a 
nearly-degenerate daughter, a result which is interesting 
in its own right. The more important implication, how- 
ever, is that if multiple decays are allowed, then the ve- 
locity factors for the rarest mode(s) to nearly-degenerate 
daughter(s) are not cancelled. 

The practical impact is that the effective branching ra- 
tio, BR e ff, for a rarer mode can be dramatically different 
than the naive BR, even by an order of magnitude. If the 
rare mode is taggable, as often happens in BSM scenarios, 
off-shell effects can alter the phenomenology significantly. 
Depending on the relative mass scales in the problem, the 
rare mode may be depleted, or enhanced to some asymp- 
totic finite value in the limit m — > M. This consequence 
is our second principal result. For the i-channel process 
we considered here, the corrections can even be negative, 
but this is not a general rule for arbitrary topologies and 
particle content. 

We don't observe massless scalars in nature, however, 
so our toy model is useful only to understand how a 
calculation performed fully off-shell can differ from the 
approximate on-shell result, simply due to other prop- 
agators which depend on q. That the decay threshold 
can experience sizeable corrections should not be sur- 
prising, as a production threshold can. Nevertheless, our 
result and its impact does not appear to be known in the 
literature. That finite-width effects can furthermore be 
enormous for rare decay modes is an important corollary, 
easily understood once one obtains the general analytical 
form of A. We note that an s-channel all-scalar process 
results in a slightly different analytical result, but iden- 
tical qualitative behavior. 

Discoverable BSM physics cannot consist purely of new 
scalar vertices, although it may involve additional inter- 
actions of a vector and two scalars (VSS) or a scalar and 
two fermions (SFF), shown in Figs.[T}: andQJl. Both types 
introduce additional momentum dependence into the ma- 
trix element. The decay matrix element for scalar decay 
to fermions (S:FF) is proportional to q 2 — (mi + m,2) 2 , 
where q 2 is the invariant mass of the resonance and rrii 
are the final-state fcrmion masses. The VSS vertex is 
proportional to the difference of the scalar particles' mo- 
menta, in turn roughly proportional to q. In general, the 
decay matrix element alters the integration over q 2 , ren- 
dering the NWA formally invalid. Our task is to see how 
much the off-shell result can differ from the naive one. 

We move on from a simple toy model to a realistic new- 
physics scenario and examine its practical phenomenol- 
ogy. For this purpose we choose the minimal supersym- 
metric Standard Model (MSSM) as our framework. Su- 
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persymmetry (SUSY) is a highly-motivated scenario, and 
there are several readily-available tools for performing 
off-shell calculations in the MSSM 7]. The MSSM spec- 
trum contains one new particle of opposite spin statis- 
tics for every particle degree of freedom in the Standard 
Model, in addition to two Higgs doublets. If supersym- 
metry exists, however, it is a broken symmetry, as any 
SUSY particles that might exist must be much more mas- 
sive than their SM partners. We might expect to find 
SUSY particles at several hundred GeV in mass, but not 
much above the TeV scale. Many SUSY scenarios, such 
as anomaly-mediated SUSY breaking (AMSB), have nat- 
ural near-degeracies in their spectrum, driven by high- 
scale vacuum expectation values for fields, which gener- 
ate the low-scale physical spectrum. Nearly-degenerate 
sparticles would still appear in cascade decays of heav- 
ier sparticles to the lightest sparticle, the dark matter 
candidate. This landscape is ripe for off-shell effects. 

One MSSM process which has no possible non-resonant 
interference diagrams is ud — > xtd (* = 1)2), with gluino 
decay to strange quark plus squark, as shown in Fig. [TJd. 
As a mixed EW-QCD process, at LHC it would be chal- 
lenging to dig out from QCD SUSY processes, unless 
SUSY were realized in a long-lived chargino scenario such 
as AMSB. It is however an excellent proxy for demon- 
strating the physics inherent in off-shell resonance ef- 
fects. That is, qg production would exhibit a very simi- 
lar matrix element effect, but the presence of QCD non- 
resonant interference would muddy the present lesson. 

Analytically at leading order in T/M and 1/s where 
allowed, the off-shell to NWA cross section deviation is: 
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with parent/daughter notation as before. In the limits 
Tfi — ► and T <C M, the first two terms go to 0, and 
the remaining correction is 0(T/M), as expected. There 
is also a residual log(s) dependence, which would not be 
noticed unless far above production threshold. 

We immediately see an important difference compared 
to the all-scalar toy model, which has no decay matrix 
element: the T/M corrections are proportional to j3~ 4 in- 
stead of (3~ 2 (differences in log terms are of minor impor- 
tance). This comes from the additional powers of q from 
the decay matrix element. This is our third principal 
result. We thus expect even more dramatic corrections 
to effective branching ratios for rare decays to nearly- 
degenerate daughters in spectra with dominant decays to 
lighter particles. These lighter particles dominate the to- 



tal width, so r to t is independent of (3 4 . It should be clear 
that the radically different behavior of off-shell BSM res- 
onances is also due to final-state masses not significantly 
smaller than the resonance mass. 

We define BR e ff as the off-shell cross section to the 
final state of interest, divided by the sum of all possible fi- 
nal state cross sections, each also calculated off-shell. For 
our example, this is three-body production of the final 
state xfssL,R- We choose the MSSM parameter space 
benchmark point SPSla 0], which also allows for gluino 
decays to lighter top quark plus stop, and bottom quark 
plus sbottom; these together generate a gluino partial 
width of 2.6 GeV, ~ 0.4% of the 600 GeV gluino mass. 
A typical width-to-mass ratio for MSSM gluinos lighter 
than about 1 TcV is 1 - 5%. 

Numerical results for BR e f / are shown in Fig. [2] for a 
scan over daughter squark mass (we assume degenerate 
1st- and 2nd-generation squarks). The shaded band de- 
lineates the region containing all 1st- and 2nd-generation 
squarks in all SPS benchmark scenarios with a heavier 
gluino. As expected, if the squarks are lighter than the 
stops and sbottoms, they dominate the total width, so 
BR e ff ss 1. However, note that for squarks much lighter 
than the gluino, the cross section receives sizeable correc- 
tions. This tree-level effect on the overall normalization 
is at a level comparable to the QCD next-to-leading order 
uncertainty. 

For sl within about 10% of the gluino mass, BR e ff 
can greatly exceed the naive NWA expectation. Where 
the BR becomes small, if it is a taggable rare mode 
(which most decays would be), then these corrections 
become important, lest we extract incorrect Lagrangian 
parameters from the relative branching ratios observed in 
data. At SPSla, this decay mode would receive almost 
a 25% correction, an order of magnitude larger than the 
T/M estimate, and greater than the residual QCD next- 
to-leading order production rate uncertainty. The NWA- 
derived ratio of this mode's branching ratio relative to 
(unaffected) decays to lighter sbottoms would disagree 
with the mass spectra measured via other methods 0]. 
Results are nearly identical for decays to s^. 

In addition to effective branching ratio corrections, we 
observe chirality selection in gluino decays. In the NWA, 
the chiral-blind Majorana-fermion gluino will decay with 
equal rates to left- and right-chiral squarks and anti- 
squarks, given equal masses (that is, up to phase space 
effects for the small electroweak mass splittings which 
typically appear). However, the initial- and final-state 
helicities are connected via the Dirac structure of the 
fermion chain. The production side of the event in our 
example involves a chargino, which selects left helicity. 
The gluino's fermion mass can flip the helicity, but does 
not give equal rates. This is our fourth principal result. 
We show it graphically in Fig. [3l 

While squark masses below a couple hundred GeV are 
already ruled out by Tevatron searches, we see that there 
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Figure 2: Ratio of effective to naive BRs (left axis) and off- 
shell to NWA cross sections (right axis) for the MSSM pro- 
cess ud — > x^ssi at the LHC (cf. Fig. (T^b)). The resonsant 
gluino has an additional partial width of 0.5,1,2,5% of its 
mass (solid, dashed, dashdotted and dotted curves) due to de- 
cays to stops and sbottoms. The first- and second-generation 
squarks lie in the shaded band for all SPS benchmark points 
with a heavier gluino. 



would be a few-percent asymmetry for nearly-degenerate 
squark masses. Left-chiral squarks prefer to decay to 
cliarginos, while 5J7(2)i-singlet right-chiral squarks can- 
not. These final states differ qualitatively and can be 
distinguished, but a detailed study is necessary to deter- 
mine what level of asymmetry would be observable. 

We have examined multiple other cases involving VSS 
and SFF vertices. Examples include sbottoms decaying 
to stop plus W boson (S:SV), and squarks decaying to 
quark plus gluino (S:FF). Their general (JofsI&nwa be- 
havior is qualitatively similar to what we find for the 
F:FS case, again depending on /3~ 4 , with minor vari- 
ations in the log terms and coefficients. Most MSSM 
particle decays would exhibit the same phenomenologi- 
cal features of our primary example. 

In summary, we investigated off-shell matrix element 
effects in scattering processes involving new heavy states 
which decay to other massive states. In general, if mul- 
tiple decay modes are allowed, the accuracy of the ef- 
fective branching ratios differs from the naive 0(T/M) 
expectations based on the NWA, even by orders of mag- 
nitude. Rarer modes may receive enormous corrections. 
If neglected, this would corrupt the extraction of model 
parameters from data. Additionally, massive Majorana 
fermions exhibit a helicity selection effect which may in- 
troduce observable asymmetries into the data. Unan- 
ticipated, such asymmetries would likely be incorrectly 
interpreted as signals of additional new physics, such as 
CP violation. Our results are based on a neglected but 
important aspect of tree-level cross section calculations. 
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Figure 3: MSSM gluino decay asymmetry to degenerate-mass 
left-chiral v. right-chiral squarks as a function of the squark to 
gluino mass ratio, in xtg production at the LHC, as discussed 
in the text. The NWA predicts exactly zero for all masses. 
The yellow (blue) curves are for multi-mode (single-mode) 
decays; line types and shaded band as in Fig. [2] 



Our results are far more generally applicable than just 
to supersymmetry. They are based on general matrix el- 
ement behavior for arbitrary renormalizable interactions, 
phase space and integration of a heavy resonance's prop- 
agator over a range in q 2 for decay to massive daugh- 
ter particles. For instance, Universal Extra Dimensions 
models are another example of cascade decays and 
very close degeneracies. 
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